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Abstract 

We consider a random trial-based telegraph process, which describes a motion on the 
real line with two constant velocities along opposite directions. At each epoch of the 
underlying counting process the new velocity is determined by the outcome of a random 
trial. Two schemes are taken into account: Bernoulli trials and classical Polya urn trials. 
We investigate the probability law of the process and the mean of the velocity of the 
moving particle. We finally discuss two cases of interest: (i) the case of Bernoulli trials 
and intertimes having exponential distributions with linear rates (in which, interestingly, 
the process exhibits a logistic stationary density with non-zero mean), and (ii) the case of 
Polya trials and intertimes having first Gamma and then exponential distributions with 
constant rates. 

Keywords Telegraph process; random intertimes; Bernoulli scheme; Polya urn model; 
logistic stationary density 
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1 Introduction 

Since the 1950s several authors investigated the (integrated) telegraph process as a realistic 
model of random motion (see Goldstein |14j , Bartlett [3J and Kac [T7j ) . Such process describes 
a motion of a particle on the real line characterized by constant speed, the direction being 
reversed at the random epochs of a Poisson process. The probability density of the particle's 
position satisfies a hyperbolic differential equation, whose probabilistic properties have been 
studied for instance by Orsingher [2U], [21] . Foong and Kanno [13], and more recently by Beghin 
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et al. [1] . See also the book of Pinsky [22] , in which many results on the telegraph process and 
its generalizations were given, and many applications were discussed. Certain one-dimensional 
generalizations of the telegraph process focus on cases in which the intertimes between two 
consecutive changes of direction are characterized by a variety of distributions. We recall 
the case of Erlang distribution (Di Crescenzo [9]), of gamma distribution (Di Crescenzo and 
Martinucci |10j). of the exponential distribution with linearly increasing rate (Di Crescenzo 
and Martinucci [II]). Moreover, investigations on the telegraph process with random velocities 
have been performed by Stadje and Zacks [21], and by De Gregorio [8], whereas a generalized 
telegraph process governed by an alternating renewal process was studied by Zacks [26]. Iacus 
[16] gave a rare example where an explicit probability law is obtained for an inhomogeneous 
telegraph process. 

In this paper we aim to study the telegraph model subject to a further source of randomness, 
by assuming that the velocity of the moving particle is driven by random trials. Precisely, we 
deal with a two-velocity random motion on the real line where, differently from the classical 
telegraph process whose positive and negative velocities are alternating, at each time epoch 
now the new velocity is determined by the outcome of a random trial. The latter follows 
one out of two schemes: the Bernoulli scheme, which acts on independence, and the classical 
Poly a scheme (cf. Polya [23j, Mahmoud [18J), where the outcome of each trial depends on the 
outcomes of the previous ones. 

The paper is organized as follows. In Section [2] we describe in detail the mathematical 
model of the motion. In Section [3] we investigate the probability law and the conditional mean 
velocity of the process. Then, in Section HI for the Bernoulli scheme we discuss the instance 
in which the random intertimes between consecutive trials are exponentially distributed with 
linearly increasing rates. In this case we obtain the probability density of the process in closed 
form, we show that it possesses a logistic stationary density, and then express the conditional 
mean velocity of the process in terms of hyper geometric functions. Finally, in Section [5j we 
show a case for the Polya scheme in which the first random intertimes in both directions 
are Gamma-distributed, whereas all remaining intertimes are exponentially distributed. We 
obtain the probability density in closed form and the conditional mean velocity as series of 
Gauss hypergeometric functions. For the reader's convenience, the paper is enriched by two 
Appendices containing definitions and formulas used in the proofs. 

2 Stochastic model 

Let {(St,Vt);t > 0} be a continuous-time stochastic process, where St and Vt denote respec- 
tively position and velocity at time t of the moving particle. The motion is characterized by 
two velocities, c and — v, with c, v > 0, the direction of the motion being specified by the sign 
of the velocity. At time Tq = the particle starts from the origin, thus Sq = 0. The initial 
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velocity Vo is determined by the outcome X\ of the first random trial. At the random time 
T\ > the particle is subject to an event, whose effect possibly changes the velocity according 
to the outcome of the second random trial. This behavior is repeated cyclically at every instant 
of a sequence of random epochs To = < T\ < T2 < T3 < • • • . We assume that the dura- 
tions of time intervals [T n , T n+ i), n = 0, 1, 2, . . ., constitute a sequence of non-negative random 
variables. Precisely, let Uk (resp. Dk) denote the random duration of the k-th time period 
during which the particle moves forward, with velocity c (resp. backward, with velocity —v). 
Furthermore, k > 1} and {D^; k > 1} are mutually independent sequences of non-negative 
and absolutely continuous independent random variables. Denoting by Z n the velocity of the 
particle during the interval [T n ,T n+ i), we assume that {Z n ;n > 0} is a sequence of random 
variables governed by the sequence {X n ;n > 1} of the outcomes of the random trials. Moreover 
we assume that the collection {Uk, Dk] k > 1} is independent of {X n ; n > 1}. 

In this paper we focus on the case in which the random trials {X n ; n > 1} are the outcomes 
of a sequence of indicator functions such that 

P{Z = c} = P{X, = 1} = J±-, P{Z = -v) = P{X 1 = 0} = (1) 

and for n > 1 



p{z n = c\g n } = p{x n+1 = i\g n } 
p{z n = -v 1 g n } = p{x n+1 = o\g n } 



b + r + An ' 



(2) 



b + r + An 

where b and r are positive constants, A is a non-negative constant, Qq = {0, Q} and Q n = 
cr(X\, . . . , X n ) for n > 1. 

We can distinguish two different cases. 

1) The case A = 0, which means that the random trials {X n ;n > 1} are independent, i.e. 
they are a Bernoulli scheme with parameter 

P=T^T- (3) 
+ r 

2) The case A > 0, which means that {X n ; n > 1} are the outcomes of a sequence of drawings 
from an urn which initially has b black balls and r red balls and then is updated according 
to the classical Polya urn scheme: the drawn ball is returned into the urn together with 
A > balls of the same color0 (X\ = 1 means that the drawn ball is black; otherwise 
we have X\ = 0.) 

In Section \5. II we discuss an extension of the above setting to the case in which {X n ; n > 
1} is the sequence of the outcomes of drawings from a randomly reinforced urn. 

The novelty of the present model of random motion is the inclusion of a sequence of random 
trials that regulates the velocity of the particle. 



We recall that, from a mathematical point of view, also in this case the parameters b, r, A of the model can 
be real. 



3 




T Tx T 2 T 3 T 4 T 5 T 6 t T \ Tj T 2 T 3 T 4 T 5 



H >+< H 




H H 

u 3 



D x D 2 



D 3 



Di 



•+« H H H 

Ui u 2 u 3 

H H 

Do 



Figure 1: Sample paths of St with (a) Vq = c, and (b) Vq = —v- 

Let Mt be the stochastic process which counts the number of epochs Tj (with i > 1) 
occurring before t, i.e. 

M t = max {i>l:Ti<t}, t > 0. (4) 

It is worthwhile to note that position and velocity of the particle at time t can thus be formally 
expressed as: 



V t = Z Mt , S t = / V s ds, t>0 



I 

Jo 



(5) 



Figure 1 shows two examples of sample paths of St , with indication of the random intertimes 
Uk and Dk, where the sequence {X n ;n > 1} takes values (a) {1,1,0,0,1,0,1,...} and (b) 
{0, 1, 1, 0, 1, 1, . . . }, respectively. 

The following stochastic equation holds: 

s n+1 =s Tk + w k , k>0, 

where k > 0} is the sequence of random variables defined by 

\ clli if Z = c 
W = < 

[ —vD\ if Zq = —v, 

and, for k > 1 and 1 < j < k + 1, 

[ cf/j if Zfc = c and Xi + A^fc_i = j — 1 
—vDj if Zfc = — u and X\ + A^^^i = k — j + 1 , 

where iVfc_i is the random variable that counts the number of random trials yielding velocity 
c among the trials going from the 2-nd to the fc-th one, i.e. 

fe-l fe-i k 

N k _x = l{z t =c} = Y; X i+i = Y. X h> k ^ 2 - (6) 



Wk 



i=l 



h=2 
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For convenience, we set also iVo = 0. In Appendix [A] we give some useful formulas regarding 
the probability distribution of N k -i conditioned on the initial velocity for the two different 
schemes. 

Throughout the paper we denote by fu k (resp. Fjj k , Fu k ) and by fjj k (resp. Fr> k , Fjy k ) the 
probability densities (resp. cumulative distribution functions, tail distribution functions) of U k 
and D k , respectively. Moreover, and fffi (resp., F^ and F^ \ F^ and F$) will denote 
the probability densities (resp., cumulative distribution functions, tail distribution functions) 
of the partial sums 

= Ux + U 2 + --- + U k , = D 1 +D 2 + --- + D k , k>l. 



3 Probability law and mean velocity 



At time t = we assume Sq = so that at time t > the particle is located in the domain 
[— vt, ct]. The probability law of St, t > 0, thus possesses a discrete component on the points 
— vt and ct, and an absolutely continuous component over (—vt, ct), which will be investigated 
in the sequel. 



Proposition 3.1 For all t > we have 

k=l 



P{S t = ct} = t—Fjj, (t) + V P{Z = c,N k = k} 
b + r * — ' 



F$\t) - F$ +1 \t) 



(7) 



where 



P{Z = c, N k = k} = { 
Similarly, we have 



P 



k+l 



b + r 



b + A 
A 



b + A + r 
A 



ifA = 
ifA>0. 



P{ St = -vt} = T^—F Dl (t) + ]TP{Z = - v , Nk = 0} F ( D k \t) - Fg +1) (t) 



where 



b + r 



P{Z = -v, N k = 0}={ 



k=l 



(8) 



(1-pf 



k+l 



r + A 



b + A + r 
A 



-i 



ifA = 
if A > 0. 



b + r \ A 

Proof. First of all, note that, for y G {— v, c}, condition St = yt implies Vq = y and V t = y 
a.s. Conditioning on Mt = k and recalling ([6]) we have 

P{S t = ct} = P{S t = ct, V t = c,V = c} 

+oo 

= Y, p {^o = c,N k = k} P{M t = k\Z = c,N k = k}. 

k=0 

Eq. ([7]) thus easily follows from (PQ) and from Eqs. (|45jl and (I50p for j = k — 1. In the same 
way Eq. (jSj) follows from (pQ) and from Eqs. (148^ and (|53p , □ 
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Let us now define the probability density of St, t > 0, conditional on initial velocity 
y € {-v,c}: 

p{x,t\y) = ^P{S t <x\V = y}, xE(-vt,ct). (9) 
We remark that, for t > 0, 

/ct 
p(x,t\y)dx = 1, y£{-v,c}. 
-vt 

The density of the particle position is 

p{x,t) = ^-V{S t <x} = -^—p{x,t\c) + -^—p{x,t \ -v), (10) 
ox o + r b + r 

where p(x, t\y), defined in Q, can be expressed as 

p(x,t\y) = f{x,t\y) + b(x,t\y). (11) 

Here / and b denote the densities of particle's position when the motion at time t is charac- 
terized by forward and backward velocity, respectively, i.e. 

f(x, t\y) = ^P{S t <x,V t = c\V = y}, 

b(x, t\y) = ^-P{S t <x,V t = -v\V = y}, 

for x in (—vt, ct). Then, for y £ {— v, c}, we can write 

f{x,t\y) = Y,fk{x,t\y), b(x,t\y) = J2 b k(x,t\y) (12) 

k=l k=l 

where, for x in (—vt,ct), 

fk(x, t\y) = —F{S t <x,Vt = c,M t = k\V = y], (13) 
b k (x, t\y) = ^-P{S t <x,V t = -v, M t = k\V = y}, (14) 

with Mt defined in (jU). Hereafter we give an expression for densities f k and b k conditioned by 
Vq = c. Note that t\c) =0. 



Theorem 3.1 For t > and —vt < x < ct, the densities Iil3\) and {1$ can be expressed 



as 



. k-2 

f k (x, t\c) = — — y~) P{N k ^ = j,Z k = c\Z = c} 
c + v ^-^ 

J=° (15) 
x f%-*- x \t - n) f f { J +1 \s - t + r*)F U]+2 (t - s)ds, k>2, 
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and 



b k (x,t\ c) = -L-{p{jV fc _i = k - 1, Z k = -v | Z = c}$\r*)F Dl (t - n) 

fe-2 

+ 1 {k>2} F i N k-i = 3, Z k = -v\Z = c} (is) 

3=0 

X /# +1) (T*) ffS^Hs - T*)F Dk _.(t - s)d S \, k > 1, 



w/iere 

. , + X , . 

T*=T*(x,t) = — ■ , (17) 

c + v 

and where P{iVfc_i = j,Zk = c\Zq = c} and P{A^_i = j, Zk = —v \Zq = c} are given 
respectively by |^5p and j/A = 0, and by |5Q|) and < f57]) i/ A > 0. 

Proof. Recalling (|13p . for i > 0, — < x < ci, and fe > 2 we have 

fk(x, 1 1 c)dx = / P{T fc G ds, Z k = c,S s + c(t - s) G dx, T fe+ i -T k > t - s, 
Jo 

< iV fc _i < k - 2 j Z = c}. 

(Case k = 1 does not give an absolutely continuous component). Conditioning on N k -i, and 
taking into account the number of time periods during which the particle moved forward and 
backward, we obtain 

fe-2 t 

f k (x, 1 1 c)dx = V / P{U {j+1) + Z^-J- 1 ) g ds, cU ij+1) - vD&'i-V + c (t - s) G dx} 
j=o Jo 

xP{U j+2 > t - s}P{N k -i = j,Z k = c\Z = c}. 

Note that conditions S s + c(t — s) = x and S s > —vs provide s > (ct — x)/(c + v) = t — r*. 
This inequality and the independence of U^ k ' and thus give 

fe-2 t 

f k (x, t\c)=J2 P{iVfc-i = j,Z k = c\Z = c} h(s, x -c(t- s))P{U j+2 >t - s}ds, 
j=o Jt ~ T * 

where h{-, •) is the joint probability density of (U^ +l "> + D^ k ~i~ l \cU^ + ^ - vD^^). Since 

*<.,,_ c(t . = (. _ *^£) /<*-;-> (-^f " 

C+V \ C+VJ \c+v 

Eq. (fT5|) thus follows recalling (fT7|) . Eq. (fT6|) can be obtained in a similar way. Indeed for 
k > 1 

6 fc (x, 1 1 c)dx = / P{T k G ds, Z fc = -v, S s - v(t - s) G dx, T fc+ i - T k > t- s, 
Jo 

< iVfe-i <k-l\Z = c}. 
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Conditioning on N k -i, and taking into account the number of time periods during which the 
particle moved forward and backward, we obtain for k > 1 

b k (x,t \c)dx = [ P{U (k) eds,cU {k) -v(t- s) edx} 
Jo 

x P{D 1 > t - s}P{N k _ l = k - 1, Z k = -v | Z = c} 

k—2 i 

+ V / P{U^ + D^-J- 1 ) G ds, cf/( J+1 ) - vD^-V _ ^ _ s ) G d:E } 
3 -=o-/o 

x P{IVj > * - s}P{iV fc _ 1 = j, Z k = -v\Z = c}. 

(Note that for k = 1 the above term PlA^-i = k — 1, Z k = —v \Zq = c} reduces to P{Zi = 
—v \Zq = c} and the sum on j is equal to zero.) For the first integral, conditions U^ k ' = s and 
c [/( fe ) _ v (f — s ) — x imply s = r*; while in the second integral, conditions S s — v{t — s) = x 
and S s < cs provide s > (x + vt)/(c + v) = r*. This inequality and the independence of 
and thus give for the above sum on j 

fc-2 . t 

P{iV fe _i = j, Z fc = - v \Z = c} h(s, x + v(t- s))P{D k _j > t - s}ds, 

3=0 Jt * 

where, as above, h(-, •) is the joint probability density of (U^ +1 > +D^ k ~ J ~ 1 \ cU^ +1 >— vD^ k ~ 3 ~^>). 
Since 

h ( a , X + V (t- s)) = -^-/S +1) ir*) fS~ j ~ 1] (« - r.) , 

we get (HBJ). □ 

Remark 3.1 Similarly to (|15p and (|16p . when the initial velocity is negative, for all i > and 
— vt < x < ct, the densities (fT3"|) and (fT4"|) are expressed as 

1 f„,„ , ( fc ) 



/ fe (x,t | - v) = ——\ P{N k ^ = 0,Z fc = c\ Z = -v}f£>(t - t*)F Vi (t* 



fc-2 



+ 1 {k>2} P{Nfc-i =k-l-j,Z k = c\Z = -v} 

3=0 

x/i i+1) (t-n) / $- j - 1 \8-t + TjF Uk _.(t-a)ds}, k>l, 

Jt-T* ' 

j fc-2 

6 fe (a;, t\ - v) = — ■ — Pj^fc-i = k - 1 - j, Z k = -v \ Z = -v} 



3=0 

t 



x fij~ j ^\r*) [ f% +1 \s - n)F Dj+2 (t - s)ds, k > 2, 



where r* is defined in (fT7|) and where P{N k _i = k — 1 — j, Z k = c \ Zq = —v} and P{N k _i = 
k — 1 — j, Z k = —v | Zq = —v} are given respectively by (I47p and (|48p if A = 0, and by ([52 
and (1531) if ^ > 0. 
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In the following proposition, conditional on positive initial velocity, the mean of process Vj 
is expressed in terms of the cumulative distribution function of T k . (Obviously, we can obtain 
a similar expression conditioning on negative initial velocity.) We remark that, by conditioning 
on the value of N k _i and using the independence of U^ +1 ^ and D^ - -? -1 ), it follows: 

fc-i 

F n\z (t \c) = J2 P i N k-i =j\Zo= c}P{^ +1 ) + D^~ x ) < t}, (18) 

3=0 

where 

P{[ /0'+i) + D (k-j-D <t}= f F { D k - j - 1} (t - s)f^ +1 \s) ds = f F^ +1 \t - s)fS- j - 1 \s) ds. 

Jo Jo 

(19) 



Proposition 3.2 For all t > we have 

+oo +oo 

Wt\Vo = c] = cF Vl (t) + c7r A ^2<l> k (t\c)+ (-«)(! -it a) ^>fc(*|c) (20) 



k=l k=l 



where ir A = and 

b + A + r 



rt 

Mt\c) = F TklZo (t\c)- I F TklZo (t-s\c)f Uk+1 (s)ds, 



fc+i 1 

t 



ip k {t | c) = F Tk \ Zo (t | c) - / F T \ Zo (t - s | c)f Dk+1 {s)ds 

Jo 



(21) 



Proof. We recall that in both cases (by independence in the Bernoulli scheme and by ex- 
changeability in Polya scheme), we have for each n > 2 

P{X n = 1 | Xt = 1} = 7TA, ?{X n = | Xx = 1} = 1 - TT A ■ (22) 

Hence, for every positive integer k > 1, recalling @ and the first of (|5|), we have 

+ 0O 

E[Vi j F = c] = E[Z Mt \Z = c} = cF Vl (t) + E l Z k ■ l{T k <t<T k+l} \Z = c]. 

k=i 

Since, for k > 1, 

E[Z k ■ l{T fc <t<T fe+1 } \Z = c] 

= cP{Z k = c \Z = c} E[l {Tfe < t<Tfc+l} \Z = c,Z k = c] 

+ (-«) p {^fc = -« | = c} E[l{ Tfc <t<T fe+1 } I = c, Zfc = -v] 

= cP{X k+1 = 1 1 X x = 1} P{T fc < t < T fc+1 | Z = c, Z fc = c} 

+ (-t;)P{X fc+1 = | X 1 = 1} P{T fc < t < T k+1 \Z = c,Z k = -v} 

/•oo 

= cir A / P(i - s < T k < 1 1 Z = c)f Uk+1 (s)ds 
Jo 

POO 

+ (-«)(1 -n A ) P(t - s < T k < 1 1 Z = c)f Dk+1 (s)ds 
Jo 

= C1T A 4>k(t | C) + (~V)(1 - TT A ) 1pk(t | C) . 
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Eq. dUD thus follows, due to (|2"Ij) . 



□ 



The following remark turns out to be useful in some particular cases. 

Remark 3.2 It is worthwhile to note that, using (|18|) and (|19|) . Eqs. (|21|) can be rewritten as 

fe-i 

0i fc (t|c)=^P{iV fc _ 1 =i|Zo=c} 

3=0 
fc-1 

c) = ^P{iV^ 1 = j\Z = c} 

3=0 

For instance, the above quantities can be easily computed when the random variables Uk 
(and Dk) are gamma distributed with the same scale parameter, since it is well-known that, 
by independence, the sums of the involved random variables are still gamma distributed. 

In the following sections we discuss some special cases arising in the two different schemes 
of Bernoulli and Polya trials, and leading to closed forms for the probability law of St- 

4 Particular case for Bernoulli scheme (A = 0) 

The classical telegraph process is characterized by exponentially distributed times separating 
consecutive velocity changes. The extension of such a model to the case of velocities driven 
by Bernoulli trials can be performed in a simple and natural way. Indeed, if we consider a 
simple telegraph process St with alternating velocities c and —v, and with alternating switching 
intensities A := (1 — p)X and ft := p/j,, then it can be proved that the marginal distributions of 
St, t > 0, and St, t > 0, are identical (where St is a telegraph process with switching intensities 
A and //, and having the same velocities of St, governed by Bernoulli trials with parameter p). 
Hence, in this case results on St can be immediately obtained by those of St- 

Therefore, aiming to discuss a non trivial case, and stimulated by previous studies (see Di 
Crescenzo and Martinucci |llj and Di Crescenzo et al. [12]) involving finite- velocity random 
motions with stochastically decreasing random intertimes, in the following we assume that the 
r.v.'s Uk and Dk have exponential distribution with linear rates Xk and fj,k. Hence, the tail 
distribution functions are 

F Uk (t) = e- Xkt , F Dk (t) = e-» kt , t > 0, (23) 

with A, fi > 0. Figure 2 shows some simulations of St in the present case, where the particle 
exhibits a kind of damped motion. Due to assumption (|23[) . f/W and D^ k \ k > 1, have 
generalized exponential densities 

f( J k \t) = k(l-e- xt ) k - 1 \e- xt , f<£\t) = k(l-e-» t ) k - 1 f ie-t lt , t>0 (24) 



F, 



(k-j-l). 



D 



(t-y) ft 1] (y)-fuu^ + u k+1 (y) 



dy, 



Fij +l) (t-y) 



Ak-j-l) 
JD 



GO - fD(*-'-V+D k+1 (y) 



dy. 
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s, 




Figure 2: Simulated sample paths of {St} in the Bernoulli scheme, exponential damped case, 
with A = /i = 1 and c = v = 1, for some choices of p. 



with corresponding cumulative distribution functions 

F^\t) = (1 -e~ xt ) k , F§\t) = (1-e-"*)* t > 0. (25) 

Hence, t/W and D( fc ) are distributed as the maximum of k i.i.d. random variables having 
exponential distributions with rates A and fi, respectively. 

By making use of the results in the previous section, in the following we obtain the proba- 
bility law of St- We start providing the discrete component. 

Proposition 4.1 Let £/& and be exponentially distributed with rates Xk and fj,k, k = 
1, 2, • • • , respectively. For all t > we have 

?{St = ct} = P ^ _, n , P{S t = -vt} - 



l-p(l-e~ Ai )' 1 J 1 - (1 -p) (1 - e-^) ' 

Proof. It immediately follows from Proposition 13.11 and Eq. (|25|) . □ 

The following theorem and corollary give the absolutely continuous component of the proba- 
bility law of St- 

Theorem 4.1 Under the assumptions of Proposition \J7T\ for all t > and —vt < x < ct, we 

have 

(l-p)^^)(e^-l) 

j{x,t\c) = k (26) 

(c + u) [pe^ + (1 - p) e ( A +A t ) T * ] 

(1 -p) \\pe XT *+^ t+T ^ + A(l - n)e AT * +2 ^ T *l 

i(i,tc=- — 2 -» ( 27 ) 

(c + v) [pe^* + (1 - p)e( x +^ T *] 



where r* is defined in (|T7|) . 



11 



Proof. Since U\ and Uk+i are exponentially distributed with parameters A and A(fe + 1), 
respectively, from ([T2"]) . (fT5l) and recalling ([23]) and ([MD we get 

+00 fc-2 



f(x,t\ c) 



C + V 



EE 

fc=2 j=0 



k - 1 
j 



p^il-pf-^ik-j -l)e 



\k-i-j 



1-e 



-£t(t-T„) 



fc-2-j 



/' A(j + l) e ^ A (^^) [l - e- A ( s -^)j ' e- A (*- s )^ +2 )d S 

Jt-T* <■ J 



M 1 -p) 

C + V 



+00 



fc-2 



e -^x;(*-i)E 



fc=2 



fe-2 



^'(1 -p) 



fc-2-j 



1 - e -^(*- T *) 



fc-2-j 



x e 



-At, 



t-T, 



X(j + l)e- A ( 4 - s ) [, 



X(t-s) L-A(t-s) _ p -Ar« 



ds 



/ip(l - g) c _^ t _ Tt) _ 2XTt , AT<I _ ^ 



c + u 

+00 fc-2 



fc=2 



j=o 



J 



E^^ET P(l-e- A V[(l-p)(l-e-^-*) 



fc-j-2 



This yields density (|26p. Eq. (|27p can be obtained from (|12p . (jT6j) and (|46p in a similar way. 
Indeed, for k = 1 we have 

6l(X)t | c) = ^i^ e -Ar,-M*-r*) 
c + u 



and, for k > 2, 

6 fe (x,t I c) 



1_ |p*-i(i _ p )jfc _ e ~ Ar *) fe 1 Ae- Ar *-^'- T *) + 
p> {I - p) k ~ j (j + 1) (l - e~ Ar *y Ae" AT * 



k-1 

j 



c + 

fc-2 

+ E 

j=0 

x / (fc - j - 1) (l - e"^ s - T *) 

Jr-t 

A 



fc-j-2 



/ie 



-fi(s-T t )-ti(k-j)(t-s) ^ g 



C + V 

k-2 

+e" A - E 
j=o 



-Ar.-M(t-T.)(l _p)jfc Ul - e - Ar *) 



fc-1 



xe 



-At(t-r«) 



k-1 

3 

t 



p(l-e" AT *) "(l-p) fc - J (j + l) 



fc-j/ 



i-i) 



e -fj,(t-s) _ e -Al(t-r») 



fc-j-2 



/^e 



^ e -AT.-A»(t-T.)( 1 _ p j 



c + t> 



fc-2 



to*- 1 + e 

j=o 



fe-1 



(j + 1)^/3 



j ok-l-j 



where we have set 

a = a(x,t) = p(l-e~ Ar *) and /3 = i) = (1 - p) (l - e"^*"^ 
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Hence, recalling that (due to binomial theorem) 



k-1 



k-2 

£ 

3=0 



from the second of (|12p we finally get 

A e - AT *-'*(*- T *)(l-p) 



(j + l)a J '/3*- 1 --' = (a + Z?)^ 1 + (k - l)a(a + /3) fe " 2 - A;a fe_1 



+00 



6(x, £ I c) 



C + V 



1 + J^(a + /5) fc - x + a Y,( k -!)(« + Z 3 ) 



fc-2 



fc=2 



fc=2 



^ e -AT.-Ai(t-T.)(l -p) 1-/3 



C + V 

which coincides with Eq. (127|) . 



[l-(a + /3)]S 



□ 



We are now able to find out the probability density of St in closed form. 



Corollary 4.1 Under the assumptions of Proposition \^T\ for all t > and —vt < x < ct, we 

have 



and 



p(x, 1 1 c) 



p(x, 1 1 — v) 



1 - p A(l - p )e( A+2 ^ r * + pe ^( t+T *) [(A + /i)e 



A-. 



/' 



C + V 



[pe^ + (1 -p) e ( A +^) r *] 2 

p e M(t+r») ^ A + ^ e Xr, (1 - p) + ^p] _ A(l - p) e ( A + 2 ^) T * 



c + u [p 6 Mt + (1 _p) e (A+M)T„] 2 

where r* is defined in (|17p . Hence, for —vt < x < ct we have 

p 1 exp{(/i-f)i-ix} 



p(x,t) 



where s := (c + v) / (X + /i). 



1-p s 



l + T ^ eX p{(/x-f)t-ix} 



2 ' 



(28) 
(29) 

(30) 



Proof. Recalling (jlip . from densities (j26[) and (j27|) we obtain ()28|) . By symmetry, density 
(f29|) can be expressed from (f28j) by replacing t* by (i — r*) and interchanging A with /i, and p 
with (1 — p). Therefore we finally get 

p(l-p){X + fi) e M*+(*+/^* 



p ( x ' ^ ( c + w ) [(i _ p ) e (A +M )T, + peM t] 2 ' 
and then Eq. ()30|) easily follows. 



< x < ct, 



□ 



Some plots of density (|30j) are given in Figures 3, 4 and 5 for various choices of the involved 
parameters. 



Remark 4.1 We can analyse the behavior of the density (130h when x tends to the endpoints 
of the state space [—vt, ct]. For t > we have 



lim p(x, t) = — — 

xl-vt 1 — p S 



1 + lV* 



2, limp(x,t) = - 

z xtct 1 — p S 



-\l 



1 + -2- e~ At 
1 ^ i-p e 
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Figure 3: Density (|30p . with t = 1, c = v = 1, A = 1, for /i = 1 (Ze/t) and fi = 2 (right), and 
p = 0.1, 0.2, 0.3, 0.4, 0.5 (from bottom to top near cc = 1). 



P(x,t) 






/ / V 20 / 




/ / /°-P 




/ / / r 10 \ 




/ / / / °7 5 






Figure 4: Same as Figure 3, with t = 10 and p = 0.1, 0.3, 0.5, 0.7, 0.9 (from left to right in 
both plots). 



The following special case follows by straightforward calculations. 

Corollary 4.2 Under the assumptions of Proposition \4^\ if Aw = fic then density 130\) be- 
comes the truncated logistic density given by 



p(x,t) 



-(x—m)/s 



s[l + e -(z-«i)/s] 

V 

1 — p 



2 ' 



-vt < x < ct, t > 0, 



m = s In 

Hence, in this case St admits of the stationary density 

e -(x-m)/s 

lim p(x, t) = — 

t— >+oo e M _j_ p—(x—m)/s 



V C 
fl A 



2 • 



x 6 



(31) 



We note that the right-hand-side of Eq. (|3ip is a logistic density with mean m and variance 
7r 2 s 2 /3. In addition, if p = 1/2 then the mean m vanishes, and the density identifies with the 
stationary p.d.f. of a damped telegraph process, as obtained in Corollary 3.3 of Di Crescenzo 
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P(x,t) 



p(x.t) 



0.6 




0.5 




0.4 




0.3 




/ \ - 2 




/ o.\ 





0.6 




0.5 




0.4 




0.3 
IS.2 




/ 0.1 





Figure 5: Same as Figure 4, for p = 0.1 (left) and p = 0.5 (right), with /x = 1, 2, 3, 4 (from left 
to right in both plots). 

and Martinucci [TT] . We remark that if A v ^ /i c then 

lim p(x,t) =0, x G K. 
t— s>+oo 

Finally, let us obtain the mean velocity conditioned by positive initial velocity. 



Proposition 4.2 Under the assumptions of Proposition for all t > we have 

+oo +oo 

E[V t |^o = c] =ce~ xt + cir A J2<Pk(t\c)+ (-v)(l-ir A ) £Vfc(<|c), 



fc=i 



fe=i 



where 



k-l 



k-l 



! (t\ c ) = xj2(j + m . yii-P?- 1 -* E E 



^=0 fii=0 



(-1) 



x [te - ^* iFi(l, 2; - A(/t + l)]i) 

-A(jfe + 1) e - A(h+1)i £T(/xi - A(/j + 1), X(k - h); i)] , 



fe-i 

^(t|c) = AE(i + i) 

i=o 



fe-1 

i 



k-j-i\fj 



■1) 



-p)*- 1 ^ E E 

fcO /i=0 

x [te-^*iFi(l,2;[^- A(/i + l)]t) 

-fi(k + 1) e~ A(/l+1)< (jd - X(h + l),/i(A; + 1) - X(h + 1); i) 
and where the function H is defined by Il5y\) in AppendixlM 

Proof. The above formulas easily follow from Theorem 13.21 Indeed, after some calculations, 
from (|18p and (|19p we find the following cumulative distribution function: 



F Tk]Zo (t\c)=x^( j+ i)( k y(i- P r^ e Er ini) 

j=o e=o h=o v - 



x ie"^ 2; - A(/i + l)]t), 
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where i-Fi(l, 2;0) = 1 and 2; z) = (e z — l)/z, for z ^ 0. Moreover, the thesis follows by 

computing the two integrals which appear in the right-hand-sides of Eqs. ()2ip , □ 



5 Particular case for Polya scheme (A > 0) 



In this section we consider a special case in which the velocity changes are governed by the Polya 
urn scheme. Let us assume the distributions of U\ and D\ are T (-^ + 1,A) andr(-| + l,/i), 
respectively, and the intertimes Ut and D^, for k > 2, are exponential with parameters A and 
(i, respectively. Therefore, U^ k ' has Gamma distribution T (-^ + k, A) and has Gamma 
distribution r + k, //), so that 



A/ 



/£>(*) 



r/A+fc ^r/A+fc-1 e ~M* 



t > 0. 



(32) 



r(i + fc) ' r(5 + fc) ' 

Notice that C/i (resp., -Di) is stochastically larger than Uk (resp., D^), k > 2. Hence, the first 
time interval along both directions is stochastically greater than the other ones in the same 
direction. This is not an unusual assumption, since in renewal theory the distribution of the 
first interarrival time is often supposed different from that of the other ones (see Chapter 2 of 
Wolff [25], for instance), as for the delayed renewal processes. Moreover, differently from the 
case treated in the previous section, U\ (resp., D\) depends on b and A (resp., r and A), which 
are the same parameters involved in the Polya urn scheme described in Section 2. 

Under the above assumptions we are able to explicitly obtain the probability law of St 
in terms of the hyper geometric function iF\(u,v; z) (see Eq. (I54D in Appendix iBl) . We first 
provide the tail distribution functions of U^ k ' and (see Eq. (157p in Appendix IBl): 



F$\t) = l 



(Xt) b / A+k e- xt 



-1*1 



l,- + k + l;Xt 



(jit) 



r/A+k e -^t 



r^ + k + i) 



iFi (l,- + k + l;fit 



t > 0, 



t > 0. 



(33) 



(34) 



In the following proposition we give the discrete component of the probability law of St- 

Proposition 5.1 Let U\ and D\ be gamma distributed with parameters (^ + 1, A) and (-^ + 1, //) , 

respectively, and let and D^, for k > 2, be exponentially distributed with parameters A and 
fj,, respectively. For all t > we have 



P{S t = ct} 



b + r~ " 1V "' ' (6 + r )r(^) 
where Fu x (t) is given by A33\) for k = 1 . Similarly, 

r(nt) r/A e-^ 



1*1 



P{S t = -vt} 



' -*W*) + 



where Fu^t) is given by (34\ ) for k = 1. 



1 b+A+r 
1, : ; fit 



A 



(35) 



(36) 
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Proof. Eq. ([35]) follows from (J7J) using (f32j) . Indeed, since 



4 fc) w - ^r i; co = / - S ) d S . 

J 



n(fc+l). 



we obtain 

P{5 t = ct} 



6 + r 
b 

b + r 



f(/,W + VP{2o = c,N k = k} \VA+k> 



b/A+k-X P -As 



1 + 

(1)* (At)" 



-A(t-) dfl 



6 + ^ + r 



6 6e" Ai (At) fe / A 



1F1 ( 1 



At - 1 



By interchanging A with fi, b with r, and £/ with D in (I35p . we immediately obtain (]36j) . □ 

For the absolutely continuous component of the probability law of St we have the following 
results. 

Theorem 5.1 Under the assumptions of Proposition HOI /or all t > and —vt < x < ct, we 

have 

C(n,t)r](n,t) 



f(x,t\c) 



t — T* 



(37) 



and 



where r* is defined in and where 

£(t*, t) := exp {— At* — /x(t — r*)} 



6 + A + r , 
iFi ( 1, j ;Ar* 



+ 



£(r*,t)?7(7-*,t) 



r/A 



(An) & M +1 [//(t-n)] 

( C + U )r(^)r(5) 



(38) 



(39) 



rj(n,t) := [At* + /x(t - r*)] 
- (An)" 1 



i*i ( 1, ^4^; An + M(t - -1- A [XT * + "<* " n)] 



.4 



6 + ^ + r 



1 ' 1 ' ~ 



AAr* 



o + A + r 



(40) 



Proof. Due to Eqs. (fl2|h (JT6J) and (f50|h and recalling ([32]h under the given assumptions we 
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have 



+00 k-2 

f(x, t\c) = — V V P{iV fc -i =j,Z k = c\Z = c} 
r. +- 7; L — ' L — ' 



C + V 



k=2 j=0 



t A fc/A+j+l (s _ t + n) fe/A+i 

-r. r(I+7+i) 



..r/A+fc-j-l/j _ \r/A+fc-j-2 



r(i + ^-i-i) 



£(n,*) 



t - T* 



+00 fc— 2 

/ b+A+r ~j " XT 



fc=2 

-OO 



fc j=0 



jfe-1 



(Ar^'Kf-n)] 



fc-i-j 



f £ THSE^ { + " r,)]*" 1 - (An)*" 1 } 



fc=2 V A Jk 

t — T* 

-(An)" 1 



I [An + M (t- r »)] 1 



1^1 I 1 



1 Fi(l,^4±I. A r.]-l 



b + A + r 
A 

AXn 
~ b + A + r 

Similarly, Eqs. (|IS]> and (j5T|) give 
1 



;Ar. + M (t-r.)J-l 6 + ^ + r 



6(x, i I c) 



P{N k -i = k-l,Z k = -v\Z = c} ^* /h F Dl (t - n 

k=l 



^b/A+k T b / A + k ~ 1 e ~XT, _ 



+00 k-2 

+ F (N k -i = j, Z k = -v\Z = c} 

k=2 j=0 



X 



^b/A+j+l T b/A+j T /A+k-j-l e -XT t -ti(t-T t ) 



r(i+i + i)r(^ + fc-i-i) 

r(An) b/A 
Ar(^)n(c + *;) 



(s - n) 



r/A+k-j-2 



ds 



r(Ar,) fc M ±2^ (1 ) fc ( An ) fc 



fc=l V A Jjfc 



+ 



k-2 



^ — 4 

/ b+A+r \ X/ 
fc=2 V A Jk j=0 



k-1 



k-l-j 



+ 



AT{ b -±£)n(c + v, 
ttn,t) ¥Z (l) fc l 

T 2^ / b+A+r \ l,| 
T * fc=2 )k K - 



e- XT *F Dl (t-r*) 



b+A+r , 
' 1 ;AT * 



'[Ar. + Mt-rOl^-CAr,)* -1 



e- Ar *F Dl (t-r*; 



+ 



£(n,i) 



[An + - r*)] 1 
(An)" 1 



f) \ A \ r 

iF! ( 1, ; An + n(t - n) - 1 



.4 



A [An + fi(t — n 
6 + A + r 



1F1 1, ; An 



A(Xn) 
b + A + r 



The proof is thus complete. 



□ 
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Corollary 5.1 Under the assumptions of Theorem I5.il for all t > and —vt < x < ct, the 
probability density of St is given by Eq. (fiflj). where 



rA(An) fc M- 1 e-^ - 
P(M|C)= (c + t,)^T(^) 



1 F 1 |l, 6 + ^ + r ;Ar,j-l 



r*(i - r*) 



p(M| -«) = — — x ^ /r + A ^ FvAr*) 



+ 



^(r*,t)^(r*,t)t 



l, & + ^ + r ;/*(*-r») 



(41) 
(42) 



r*(t - r*) 

where r*, £ and are defined respectively in Eqs. (77] ), < fgg|) and |^0[ ), and where 



£(t*, i) := exp {-At* - - r*)} 



(c + ,)r(^)r(A) 



r?(r*, i) := [At* + - r*)] 
-[^(t-n)]" 1 



( i, ^±4^; An + M (* - ^ - 1 - A [Ar * + * - T * )] 



A 



b + A + r 



A 



b + A + r 



Proof. Eq. (|4ip immediately follows from (jlip . From Remark 13.11 by interchanging \i with 
A, (t — r#) with r#, & with r, and D with J7 in o(x,i | c), given in (f38l) . we get the following 
density: 



+ 



b+A+r , 
1 F 1 ( 1, m(t-n) -1 



t — r* 



Similarly we obtain the density b(x,t \ — v) by interchanging \x with A, (t — r*) with r*, 6 with 
r and D with ?7 in /(x, t\c), given in (f37|) . so that 



b(x, 1 1 — v) 



Hence, Eq. (02} follows by (fTTI) . 



□ 



Various plots of density p(x,t) given in Corollary 15. II are shown by Figures 6, 7 and 8. 



Remark 5.1 From the formulas given in Theorem 15.11 we have 



g(n,t)n(n,t) = 
t*-s>o t — r* 



lim gfr*'*)^*'*) = 0, 



Um g(r>,t)^,t) =0> lim g(n,t)n(r„t) =a 
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p(x,t) 











0.04 


0.02 





1.0 -1.0 




Figure 6: Density p(x,t) for Polya scheme with densities (I32p . with t = 1, c = v = 1, 
A = /i = 1, A = 2, for 6=1 (Ze/t) and 6 = 2 (right), and r = 1, 2, 3, 4 (from bottom to top 
near x = 0). 



p(x.t) 



U.UD 
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P(x,t) 


0.06 
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0.04 




— — 0.03 




0.02 




0.01 





1.0 -1.0 



Figure 7: Same as Figure 6, with r = 1 and A = 0.4, 0.6, 0.8, 1 (from bottom to top near 
x = 0). 

Hence, some calculations allow to study the behaviour of the density p(x,t) in proximity of 
the endpoints —vt, ct: 

6r^(/ii) r / A - 1 e-^ 



lim p(x, t) = ■„ 

x-\-vt FK ' ; (b + r)(c + v)AT(r/A + 1) 



and 



limp(a;, t) 

x^ct 



brX(Xt) b / A ~\ 



-\t 



+ r)(c + v)AT(b/A + 1) 



6 + ,4 + r 

i, — ^ — ;m* 



1^(1,^4^^1-1 



.4 



Note that in this case a stationary distribution does not exist. Indeed, 

lim P{S t = ct} = 0, lim P{S t = -vt} = 0, lim p(x, t) = 0, 

t— S-+0O t— > + 0O t— > + 00 

since (see Equation 13.1.4 of [T|), as \z\ — > +oo, 

iJfi(o,6;«) = ^eV^l + Ond- 1 )] for Re(z) > 0. 
r(a) 

Finally, the mean velocity subordinated to the positive initial velocity can be expressed in 
the following form. 
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p(x,t) 
0.07 r 




10 -10 




Figure 8: Same as Figure 6 (from bottom to top near x 



with t = 10. 



Proposition 5.2 Under the assumptions of Theorem \5.1l for all t > we have 



Wt \V = c} = cFu.it) + ctt a <t>k(t I c) + (-«)(! - VT A ) I c) 



(43) 



fc=i 



fe=i 



with 



k-l 



</>k(t\c) = J2~P{Nk-i=3\Zo = c} 

3=0 

x [G(r/A + fc-j-l,/x,6/^+j + l,A;t)-G(r/^ + A ; -i-l,/x,6/A + j + 2,A;t)] , 
fe-i 

^fe(i|c) = ^P{iV fe _i=j|Zo = C } 
i=o 

x [G(6/4 + j + 1, A, r/A + k-j- *) - G(6/,4 + j + 1, A, r/4 + fc - j, /x; t)] , 

(44) 

where P{N^_i = j \Zq = c} zs given by f^ffi ) and i/ie function G is defined in Appendix W[ by 
Eq. (S3). 



Proof. Eqs. (I43p and (I44p easily follow from Proposition [321 Remark 13. 21 and Eq. (|58h since, 
for k > 1, C/fc+i an d -Dfc+i are exponential with parameter A and /z, respectively, and and 
are Gamma-distributed with parameters (b/A + fc, A) and (r/A + A;, /i) respectively. □ 



5.1 A suitable extension 

In this section we propose a suitable extension of the model based on Polya scheme by replacing 
the classical Polya urn with the following randomly reinforced urn (cf. Aletti et al. [2], Berti 
et al. [5], [BJ, Crimaldi [7]). An urn contains b > black balls and r > red balls. At each 
epoch n > 1, a ball is drawn and then replaced together with a random number of balls of the 
same color. Say that B n > black balls or R n > red balls are added to the urn according 
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to whether X n = 1 or X n = 0, where X n is the indicator function of the event {black ball at 
time n}. We assume, for each n > 1, 

(B n , R n ) independent of [X x , B u R u X 2 , B n „ u R n _i,X n ) . 

According to this new model, formula ([2]) must be replaced by 

b + Sfc=i BkXk 



p{z n = c\g n } = p{x n+1 = i\g n } 
p{z n = -v | g n } = p{x n+1 = o\g n } 



b + r + Efe=i (B k X k + R k (l - X k )) ' 

r + ELiRk(i-x k ) 



b + r + ELl i B kXk + RkO- ~ X k )) ' 
where 

g = {0, S7}, g n = a[Xi,Bx,Ri, . . . ,X n , B n , R n ). 

It has been proved that, if B n = R n for each n (cf. Berti et al. [6] and Crimaldi [7]) or, more 
generally, i(j B n = R n for each n (cf. Berti et al. [5]) P- 538), then {X n ) n >\ is a sequence of 
conditionally identically distributed random variables (which is generally not exchangeable) 
and so we have, for n > 2, 

b + B! 



P{X n = 1 1 X x = 1} = E 



6 + Bi + r 



P{X„ = | Xl = 1} = E 



b + Bi+r 



Using the above equalities, instead of (|22p . we can obtain a formula for the mean velocity similar 
to pop . However, in this case the computation of the conditional probability distribution 
of iVfc_i needed in Eq. (|18p is not an easy task, and possibly will be the object of future 
investigations. 

A The conditional probability distribution of N^^i 
A.l Bernoulli scheme (A = 0) 

Using the independence of {X n ;n > 1} and setting p as in ([3]), we have for k > 1 and 
< j < k - 1 



PjiVfe.! = j \Z Q = c} = P <\J2 X h = 3 \X\ = 1 j = ^ ^. -P) 



k-i-j 



and so 



P{iY fc _ 1= j, Z fe = c|Z = C }= ( k 1 V' +1 (l-p) fc - 1 - J ', (45) 
which simply becomes P{Zi = c | Zo = c} = p for k = 1 and j = 0. Similarly, we have 

P{iV fe _! = j, Z fc = -« | Z = c} = f k ~ ^ (1 - p) k ~j . (46) 



2 d 

The notation = denotes equality in distribution 
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Conditioning on Zq = —v, we get 

P{iVfc_i = k - I - j, Z k = c\Z 

and 



-v] 



k-l 



P{N, 



k-l 



k - 1 - j, Z k 



-V I Zq 



-V} 



3 

k-l 



p k - j (l-p) j 

fc - 1 -j(i- p y'+ 1 



(47) 
(48) 



A. 2 Polya scheme (A > 0) 

Using the exchangeability of {X n ; n > 1}, we have for k > 1 and < j < k — 1 

P{JVt _ 1= ,| Zo = c} _,'*-A r(^) r( J + ^)r(*-i-; + 4 ) 



j 7r(^)r(j) r(t-i + !±4±i) 

fc-lAA + AA /r\ Wb+A+r 
3 ) \ ^~)\A~)k- 3 -i 



A 



k-l 



(49) 



where (a)o = 1 and (a)j = a(a + 1) • . . . ■ (a + j — 1) = T(a + j)/T(a) (the ascending factorial). 
For k > 2, the above formula corresponds to the quantity P{Ylh=2 = J I X\ = 1} 5 i-e. to the 
probability of obtaining j black balls in k — 1 drawings in a Polya urn scheme starting from r 
red balls and b + A black balls. Hence we get 

P{N k -i =j,Z k = c\Z = c} = P{X k+1 = 1 1 N k _ t = j, X 1 = 1} • P{JV fc _i = j | Z = c} 

-P{7Vfc_l = J | Z Q = C\ 



b + r + Ak 

k-i\T{ b -±j^)r(j + i + b -±£)r{k-i-j + %) 
3 



■b+A 



^;r(i)r(/ c + — 



k-l\ (b + A 

3 



A 



1 \AJk-j-i 



b±A+r\ 
A ) 

fb + A + r 
V A 



(50) 



which simply becomes ~P{Zi = c | Zq = c} = gr^4^ for = 1 and j = 0. Similarly, we have 
P{iV fe _ x = j, Z fc = -« | Zq = c} = P{X fc+1 = I N h - X = j, X 1 = 1} ■ P{JV fc _i = j | Z = c} 

*-i y(^)r(j + ^)r(fc-j4) 
i >/ r(^)r(^)r(^ + ^±r) 

6 + ^ + r 



A; — 1\ /& + A\ rr 
3 J \~A~Jj\AJ, 



A 



(51) 



Conditioning on Zq = —v, we get 
P{N k ^ = k - 1 - j, Z k = C | Zq 



-v} 



k-l\ (r + A 
3 



A 



k-j 



fb + A + r 
V A 



(52) 
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and 

P{iY fc _x = k-l-j, Z k 



-v I Z n 



-v} 



k-l\ f r + A 
j 



A 



b 

A) h 



k-j-i 



b + A + r 
A 



(53) 



B Auxiliary functions 

For the reader's convenience, now we provide some formulas used in the proofs. 

The confluent hypergeometric function (or Kummer's function) used throughout the paper 
is so defined (see [I] for details) 

{u) k z k 



iFi(u,v;z) = ^2 



k=0 



(v) k k\ ' 



(54) 



We start recalling the relation 

i-Fi(a, b; z) = e z iFi(b — a, b; — z) or, equivalently, iF\(a, b; — z) = e~ z \F\{b — a, b; z) . (55) 
An useful integration formula is Equation 3.383.1 of |15j : 

f V u ~\t - yy-^&y = T ^ )T } V \ t^ v -\F x {u, » + u; fit) (56) 

Jo r(/x + v) 

when Re(fi) > and Re(z^) > 0. 

Let X and Y be independent gamma-distributed random variables with parameters (a, ji) 
and (/?, A), respectively. By ()55[) and the integration formula (|56p . we find 

(MY 



Fx(t) 

T{a + 1) 

and so, for a, P, fj,, A > and t > 0, 



e~"* 1*1(1,0 + 1;^), H>0, t>0 



(57) 



G(a,n,(3,\;t) :=P(X + Y<t) = [ F x (t - y)f Y (y)dy 

Jo 



+oo 



(/it)' 



(58) 



h=0 



T(a + P + h + l) 



1 F 1 {p,a + P + h + l;{pi-X)t). 



Note that an alternative expression of G was given in Moschopoulos [19] as series of suitable 
integrals. Analogously to (J3HJ) , by (JoT)]) we can express the function 



H{a,p-t) := [ (t-y)iF 1 (l,2]a(t-y))e- 
Jo 

in terms of the hypergeometric function as 



e~ py dy, a, PeR, t>0 (59) 



r t 



- [ 1 F 1 (l,2;-Pt)-e- at 1 F 1 (l,2;(a-P)t)] , a ± 0, P G 



a 

t r 

p 



e-P- iFi(l,2;-0i) 



I 2 



a = 0, /3 / 



a = 0, /3 = 0, 



where we recall that 2; z) reduces to 1 when z = and to (e z — l)/z when z ^ 0. 
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